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ABSTRACT 

This paper studies the motion of a particle whose tune is inside and 
near a linear half-integer stopband. Results are found for the tune and 
beta functions in the stable region close to an edge of the stopband. It is 
shown that the eigenvalues and the eigenfunctions of the transfer matrix 
are real inside the stopband. All the results found are also valid for small 
accelerators where the large accelerator approximation is not used. 



1. Introduction 

Inside a linear-half integer stopband the particle motion can be unstable and grow expo- 
nentially. A result is found for the growth parameter g which is the imaginary part of the 
tune. It is shown that the real part of the tune inside the stopband is constant at the 
value q/2, where v = q/2 is the center of the stopband, q being some integer. The proof 
given does not depend on perturbation theory, and the result follows from the symplectic 
properties. The eigenvalues and eigenfunctions of the transfer matrix are shown to be real 
inside the stopband. In the stable region near an edge of the stopband, the tune varies 
rapidly and the beta function becomes infinite as the unperturbed tune approaches the 
edge of the stopband. Results are found for the tune and beta function in the stable region 
near an edge of the stopband. It is found that the beta function becomes infinite inversely 
as the square root of the distance of the unperturbed tune from the edge of the stopband. 
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The basic equations used are also valid for small accelerators, where the large accelerator 
approximations is not used, and all the results found are valid for small accelerators. 



2. Results when the tune is not in the stopband 

Before treating the interesting case where the tune is inside the stopband, it will be helpful 
to first treat the case where the tune is not inside the stopband. It will then become clear 
where the perturbation solution breaks down, when the tune is inside the stopband, and 
how the perturbation solution can be repaired. 

It will be assumed that in the absence of the perturbing fields, the tune of the particle 
is i/Q and that the motion is stable when vq is close to q/2, where q is an integer. 

It is assumed that a perturbing field is present which is given on the median plane by 

ABy = -G{s)x (2-1) 

G{s) is periodic in s and contains the field harmonics that can excite the stopband around 
uq = q/2. 

Introducing 77 defined by 

V = x/(3'/\ (2-2) 

where P is the beta function of the unperturbed field, the equation of motion can be written 
as 

f = 413^1^ ABy/ Bp ^2_3) 
/ = -vlp^Grj/Bp 
Bp = pc/e, d6 = ds/vQf3 

Eqs. (2-3) are valid for large accelerators, and a small change is required^ to make them 
valid for small accelerators (see section 7). However, the final results found below are valid 
for small accelerators that require the use of the exact linearized equations. This is shown 
in section 7. 
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Eq. (2-3) can be written as 



. (2-4) 

Because b{9) is periodic a solution for rj will have the form 

ri = exp (ii^sO) h (9) (2-5) 

where h{9) is periodic. It is assumed that the tune uq will change to Ug because of the 
perturbing field. Thus 77 can be assumed to have the form 

rj = As exp (ii's^) + exp {ivj.9) 

r+s (2-6) 

Vr = Vs-\- n 

where n is some integer but n ^ 0. For a zero perturbing field the solution for rj is 
r] = Aex.p{ii'Q9). Thus for small perturbing fields it can be assumed that 



Ar As for r ^ s 

Putting Eq. (2-6) into Eq. (2-4), one obtains a set of equations for the Aj. 
{yl — i/q) Aj. = 2vq ^ brfAf 



(2-7) 



r 

1 /-^^ (2-8) 
i>rr = d9 b{9) exp {—ii'r9 + iuf9) 



27r _ „ 

IVj. = I/g + 77, 

Eqs. (2-8) are a set of homogenious equations for the Ay. The condition required for 
a solution to exist is that the determinant of the coefficients of the A^. should vanish. This 
condition will determine Vg, and then Eq. (2-8) can be solved for the A^ in terms of Ag. It 
is more convenient to solve Eqs. (2-8) using an iterative perturbation procedure. For the 
initial guess for r] in this iterative procedure, one can assume 

7] = As exp {ivs9) (2-9) 

One can put this result for -q in the right hand side of Eq. (2-8), and solve for Ay. which 
gives 

iv^ -vl)Ar = 2vQbrsAs 

(2-10) 

l/^ = Vs + Tl 
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For r = s, Eq. (2-10) gives 



which determines u^,. One finds 



which can be written as 



{^s - ^o) = 2t'o6ss^s 

i's-J^o = bss 
i^s-J^o = bo 



(2-11) 
(2-12a) 



bo 



bo 



1 

2^ 
1 

4tt 



2n 



dOb (d) 



/3G/Bp 



(2.12b) 



This is the well known result^'^ for the first order tune shift due to a gradient G. The 
solution for rj to first order can be found using Eq. (2-6) and Eq. (2-8) 



T] = As I exp {iiysO) + ^ ^ brs exp {iurO) 



(2-13a) 



Ur = l^s + n, 
which can be also written as 



n 7^ 



ij = As exp (wsO) < 1 + 



2z/o6n 



n{n+ 2^0) 



exp {in9) 



(2-13b) 



1 

bn = — / dsPGexp {—inO) /Bp 
47r Jo 

One can see from Eq. (2-13b) that the above solution will not be valid when vo is close to 
a half-integer for the denominator (n -|- 2uo) can become very small. 



3. Results when the tune is inside the stopband 

It is assumed that vo is near vo = q'/2, q being an integer. The stopband is defined as the 

range of vo for which the tune, Vsi in the presence of the perturbation given by Eq. (2-2) 
has a non-zero imaginary part. One can write Vs as 

i^s = i^sR - ig (3-1) 

It will be shown that inside the stopband, where g ^ 0, then 

i^sR = q/2 (3-2) 
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This may be shown as follows. Let be the phase shift for a period, and = I'nvg where 
the period has been assumed to be one turn or 27r in ^. Let T be the transfer matrix for 
one period. Then one has 

cos// = ^(Tn+T22), (3-3) 
and one sees that cos yt, is real even inside the stopband since the Ty are real. One also has 

cos yt, = cos 27r {ugR — ig) = cos {27ri'sR) cosh {27rg) + i sin {2T^Vgji) sinh 2'Kg (3-4) 

In order for cos ji to be real, one has to have either ^ = or, if 0, 2ttVsR = "^^t^i where 
n is some integer. Thus, inside the stopband where g ^ 0, 

i^sR = n/2 (3-5) 

In order to have continuity when the perturbation goes to zero, one has n — q and 

UsR = q/2 (3-6) 

For the range of vq which is inside the stopband, VgR remains constant at VgR = q/2 while 
g goes, as will be seen below, from g = at the edge of the stopband to a maximum value 
in the middle of the stopband, where uq — q/2. One may note that this result does not 
depend on the use of perturbation theory and is valid even for large perturbations. A proof 
of the result, based on perturbation theory, is given below. 

Now let us return to the problem of computing the growth factor, g, using perturbation 
theory. The iterative perturbation approach used in section 2 appears to breakdown when 
vq is close to q/2. One sees from Eq. (2-10) that when i^q ~ q/2, then one of the Ar 
becomes comparable to As, and this is the Ar for which Ur — — T^r — ~<?/2. Thus 
in the above iterative procedure for finding r/, one will assume for the initial guess for -q, 
instead of Eq. (2-9), 

■q = As exp {ivsO) + Ag-exp {iv^) 

(3-7) 

Us = vs-q 

Then Eq. (2-10) is replaced by 

ivl - Ug) Ar = 2uobrsAs + 2uobrsAs 

(3-8) 

I'r = + IT' or i'r = ^ + n = i's — q + n 
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For r = s and r = s one obtains 2 equations for Ag and Aj 

{^s - ^o) = '^^obssAg 

{4 - i^i) Aj = 2i^oh,sAs (3-9) 

In Eq. (3-9), it has been assumed, for simplicities sake, that bgs = 0. This can be 
accomphshed by redefining uq to be + Bq. 

In order for Eqs. (3-9) to have a solution, one must have 

i^s-^o) (^l-^o) =4i^olM' 
i/g = Us-q 

Eq. (3-10) determines Vg- If one writes Us — i^sR ~ one finds 



(3-10) 



(j^sR - igf - J^o] [{j^sR-Q-igf =4t'o|Azyp (3-lla) 



where Au = 6, 



Aiy = — [ ds Pexp{-iqe)G/Bp (3-llb) 
47r Jo 

\Au\ will turn out to be the half width of the stopband. Assuming that vq and VsR ^ire 
close to g/2, then Eq. (3-lla) can be rewritten as 

(j^sR -ig- J^o) (q - ^sR + ig- t^o) = I Ai/p (3-12) 

The imaginary part of the left hand side gives the equation, 

g{2usR-q)^0 (3-13) 
Inside the stopband where g ^0, one gets 

UsR = g/2 (3-14) 

Using Eq. (3-14) for UgR, the real part of Eq. (3-12) gives 

{q/2-uof+g^ = \Auf 

r,. ,2 . . x2lV2 (3-15) 

Eq. (3-15) shows that the growth parameter g has a maximum at uq = q/2 where g = \Ai/\, 
and decreases to zero at uq = g/2 ± lAu]. The stopband width is then 2|Aiy| and extends 
from g/2 — |Ai/| to g/2 -|- lAu]. Thus |Ai/| is the half-width of the stopband. 
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Now let us find the solutions for 77 that will give the particle motion inside the stopband. 
To lowest order, 77 is given by Eq. (3-7) and inside the stopband 

rj — As exp (ivsO) + Aj-exp {ivg6) 

Us = q/2 - ig (3-16) 
Us = -9/2 - ig 
As and A-g are related by Eqs. (3-9) which gives 

Jig — 2 

{q/2 + ig) -vl 
{q/2 - z/o) + ig 

where it has been assumed that uq is close to q/2. Eq. (3-17) can be written as 

As = As exp [-i {61 + 62)] 

61 = phase Au (3-18) 

62 = phase [{q/2 - uq) + ig] 

where bgs = |Ai/| exp(— z^i), and lAu] = g"^ -\- {q/2 — vof were used. 
Putting these results for Ag into Eq. (3-16) one gets for r] and x 

r] = As exp {gO) cos {qe/2 - {61 + ^2) /2) 
6\ = phase Av 

(3-19) 

62 = phase [{q/2 - i^q) + ig] 

where the multiplying constant exp[—i{6 + 62) /2] was dropped. Eq. (3-19) give to lowest 
order the two solutions inside the stopband corresponding to whether g is positive or 
negative. The first order correction to 77 could also be found using Eq. (3-8) and the result 
is 



7] = As exp {gO) 



cos {qe/2 - {61 + 62) /2) 



(3-20) 
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The results found above for the stopband width, the growth rate and particle motion in 
the stopbands will apply equally well to small accelerators where the exact linear equations 
have to be used, without the approximations used for large accelerators, as shown in section 
7. 



4. The realness of the eigenvalues and eigenvectors inside the stopband 

It will be shown in this section that the eigenvalues and the eigenvectors of the one period 
transfer matrix are real when vq is inside the stopband. If T(s) = T(s + L, s) is the one 
period transfer matrix, then the eigenvalues and eigenvectors satisfy the equation 

t (s) x{s) = Xx (s) 

(4-1) 



X 



X 
Px 



The symbol x is used to indicate both the column vector x and the first element of this 
vector. The meaning of x should be clear from the context. The eigenvalue A is related to 
zv- value by 

A = exp (i/j,) 

(4-2) 

= 27TI' 

where, for simplicity, the period is assumed to be one turn or 27r in ^. 

In section 3 it was shown that when z/q is inside the stopband, v can be written as 
V = q/2 — ig one finds for x, when is inside the stopband, 

A = exp {iqiT + 2TTg) 

(4-3) 

A = (-l)^exp (27r^) 

and the eigenvalue. A, is real when vq is inside the stopband. 

Now let us consider the eigenvectors of T. If A, x are an eigenvalue and eigenvector of 
T, it is known that lA is an eigenvalue of T, as T is symplectic. Since T is real, A*, x* are 
also an eigenvalue and eigenvector respectively of T. It appears that one has 3 eigenvalues, 
A, 1/A and A* instead of just two. Two of the 3 apparent eigenvalues must be equal. 

When vq is not in the stopband then Vs is real, 1/A = exp(— z/x) = A*, and the second 
solution is A*, a;*. 
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When vq is in the stopband, then v has an imaginary part, v = q/2 — ig. Then 
1/A 7^ A*. To avoid having 3 eigenvalues, one needs to have A* = A and x* = x, or X and 
X are real when vq is in the stopband. 

One should note that the statement x is real means here that x can be made real by 
multiplying by a constant. This is because an eigenvector can be multiplied by a constant 
and still be an eigenvector of T. Thus the eigenvector may not appear real because of the 
presence of a constant complex multiplier, but can be made real by multiplying it with a 
constant. 

In section 3, where rj was found by an iterative perturbation procedure for z/q in the 
stopband, the solution found, Eq. (3-17), is real as required by the above theorm, whereas 
when uq is not inside the stopband, the solution for rj found Eq. (2-13) is complex. Since 
the solution for rj was specified to have the form given by Eq. (2-6), then x = rj/P^^'^ is 
the first element of the eigenvector, x, of T. 

5. Tune near the edge of a stopband 

In this section, a result will be found for the tunc in the stable region outside the stopband 
but close to one of the edges of the stopband. It will be shown that close to the edge of a 
stopband, 

\i^ - q/2\ = {2\Au\ \uo - (5-1) 

v is the tune in the presence of the gradient perturbation, Vg is the edge of the stopband, 
Ve = 9/2 ± lAi/j. |Ai/| is the half-width of the stopband. Eq. (5-1) shows that when uq is 
close to an edge of the stopband, Vg, v varies rapidly with vq, and the slope of the v vs. 
vq curve is vertical at = u,,. 

To find p in the stable region outside the stopband, where pq — q/2\ > |Az/|, one goes 
back to the derivation given in section 3 for v inside the stopband, starting with Eq. (3-12). 
Eq. (3-15) shows that for ji/Q — q/2\ > |Ai/| the only acceptable solution is ^ = 0, and Eq. 
(3-12) can be written as 

{1, - uo) {\u - q\ - po) ^ \Ap\'^ (5-2) 

where we have put = p- 



10 



Assuming that vq is just below the stopband edge = q/2 — |Ai/|, put vq = i^e — ^ 
and V = q/2 — 5 into Eq. (5-2), where e and 5 both approach zero as approaches the 
stopband edge. We find 

5 = {e(e + 2|Az/|)}i/2 ^5_3) 

The top edge of the stopband can be treated in the same way and both results can be 
combined into the one result 

\v-q/2\ = {\vQ-v^\{\vQ-Ve\ + 2\I^v\)fl'^ 

(5-4) 

Ve = q/2 ± 

Very close to the stopband edge, — ^ |Az^|, one finds 

\v-q/2\ = {2\/^v\\vQ-Ve\fl'' (5-5) 
Thus, as z/Q approaches a stopband edge, v approaches g/2, and dv/duQ become infinite 



6. The beta function near the edge of the stopband 

In this section, a result will be found for the beta function in the stable region outside the 
stopband, but close to one of the edges of the stopband. It will be shown that close to 
edge of a stopband 

[(/? - /5o) //3o]^ax = [2| AH/I^o - ^e|]'/' (6-1) 
Ve is the edge of the stopband, = q'/2± lAi/]. lAi^l is the half-width of the stopband. vq^ 
(3o are the unperturbed tune and beta function. Eq. (6-1) shows that when uq approaches 
the edge of the stopband, {(3 — (3q becomes infinite like l/|z^o ~ ^e\^^'^- 

The beta function (3 can be found from the solution for the rj function which has the 
form given by Eq. (2-5). It will be shown below that (3 is given by 



/3 


V 


^2,/ 1 
\ ?7 2 




1 

^ 2^ 


1-2-K 




/ de 

Jo 



1 \ 1 r" .„ 1 ^^'^^ 



z/, (3 and uq, (3q arc the perturbed and unperturbed values of the tune and the beta function. 
Thus P can be determined from the rj function, if the perturbed tune is known, and in this 
case the tune was found in section 5. 
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To derive Eq. (6-2) one notes that 



(6-3) 



and thus 



f3 = Dl3o \r)\^ 

where D is some constant. The constant D can be determined using the relationship 

---TdO^-^ (6-4) 

which when written in this form, wiU hold also for small accelerators^ when the perturba- 
tion does not change the closed orbit (see section 7). This gives 



{/3//3o)^{i^o/i^)\v\Hl/\v 



(6-5) 



One now proceeds to find the 77 function in the stable region near an edge of the stopband. 
To lowest order, 77 is given by Eq. (3-7) as 

T] = As exp (ivsd) + ^jexp (iugO) 
Ag is given by Eq. (3-9) which can be written as 



(6-6) 



A^ 



|Az/ 


exp 
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A.. 



(6-7) 



5i = ph (Azy) 

Assuming vq is below the edge of the stopband one puts Vs = q/2 — S and vq = Vg — e, 
Ue = q/2 - \ Aiy\. Then 



|A;y[exp(-/di) ^ 
d -I- e -I- lAi/] 

The r] functions is then given to lowest order by 



(6-8) 



1 + 





Av 




\Au 


+ e + 5 



exp \% {-qd - Si)] 



T] = As exp (ivsO) 

e=\uQ-Ue\, Ue = q/2±\Au\ 

5=\u- q/2\ = ||j^o - l^e\ (ko - l^e\ + 2\Au\) } 



(6-9) 



1/2 
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Eq. (6-9) has been written so that it holds both above and below the stopband, and Eq. 
(5-4) was used to replace 5. 
One then finds 



|r/p = \As\^ (1 + 2Ccos {q9 + 6i) + C^) 
C = 





\Aiy 




\Au\ 


+ 






+ 


\u-q/2\ 



1 1 /■^'^ 1 



\ri\^/ \As\^ 27rjQ 1 + 2C cos (qO + Si) + 
1 1 



(6-10) 



(6-11) 



|A,|2l-C2 

One can now use Eq. (6-2) to find P 

£ = ^ (1 + C2 + 2Ccos {qe + Si)) (6-12) 

One may note that when i^q is at the edge of a stopband, ~ = 0, u — q/2 = and 
C = 1. Since one is interested here in the rather large effects due to the 1/(1 — C^) factor 
near the edge of the stopband, one can put i^o/iy — 1 with only a small error. One finds 
for P 

2C {C + cos {q9 + Si)) 

[P - PO) /PO = ^ _ ^2 

2C 

{{P-Po)/Po)^^^Y:-^ (6-13) 
((/?-/3o)//3o)^,,= 



2|Az> 






+ I 


v-q/2\ 



using the result \i' — q/2\ = {{i/q — i'e){\i^o — i^e\ + 2| Ai/I)}-*^/^ one gets 

HP - Po) /Po)ra^ = —J, '^-^P ,.-,1/2 (6-14) 

Very close to the stopband, \vq — f^, <^ oii® finds 



2|Az/| 



1/2 



{{P - Po) /Po)ras. = r^-^ (6-15) 

Eq. (6-15) shows that as vq approaches the edge of a stopband {{P — /9o)//3o)max becomes 
infinite like {|z^o ~ ^e\}~'^^'^- 

Table 1 below shows how the distortion in the beta function depends on how far uq is 
from the edge of stopband. The table shows that to keep {{P — Po)/Po)max less than 10%, 
then uq has to be about lOAu away from the stopband edge. 



13 



Table 1: ((/? — Po)/Po)max versus ji^o — zvel/l^^^l as computed from Eq. (6-14). 



\uo - Ue\/\Aiy\ 


.5 


1 


2 


4 


6 


8 


10 


((/3-/3o)//3o)max 


1.24 


.73 


.41 


.22 


.15 


.12 


.095 



7. Comments on the small accelerator results 

All the final results found in this paper will also hold for small accelerators where the exact 
equations of motions have to be used. The exact linear equations have the form^ 

dx 

— = Anx + AuQx 

as 

A ^ A (7-1) 

— = A21X + A22qx ^ ' 

as 

Qx = Px/P 

In the large accelerator approximation, it is assumed that An = A22 = 0, A12 = 1, and 
Qx — dx/ds. The coefficients Aij are given in reference 1. In particular, A12 is given by 

A,, = Jiii/^^. (7-2) 

where the right hand side is evaluated on the closed orbit. 

Although for small accelerators, the differential equations for x and /3 are different from 
those usually used for large accelerators, it has been found^ that the linearized differential 
equation for -q = xjP^I'^ is not much different. One finds 

^ uiPiABy (7-3) 
A12 Bp 
Bp = pc/e 

with the difference that 6 is now defined by^ 

d9 = Ai2^ (7-4) 

^OPO 

where i/q, Pq are the unperturbed tune and beta function. The A12 is evaluated on the 
unperturbed closed orbit. 
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The relationshp between v and (3 is somewhat different for small accelerators, and is 
given by-*^ 



1 



where C is the circumference of the accelerator. Using Eq. (7-4) this can be written as 

- = ^r^^ide (7-6) 

where A12 is the A12 coefficient evaluated on the unperturbed closed orbit. If the per- 
turbation being considered does not change the closed orbit, then A12 = A12 and one 
has 

- = -r^dO (7-7) 

Eq. (7-7) now holds for both large and small accelerators, provided the perturbation does 
not change the closed orbit. 

The results found in this paper are based on Eq. (7-3) for r] and Eq. (7-7) that relates 
u and /?. Although, the perturbation term in Eq. (7-3) now has the extra factor I/A12 for 
small accelerators, this factor of I/A12 disappears in the final result when one goes from 
the variable 9 to the variable s according to Eq. (7-4). Thus the result for the stopband 
width, Au, as given by Eq. (3-llb) is valid for both large and small accelerators. Keeping 
the above equations in mind, one can go through the derivations on the previous pages, 
and show that the final results are valid for both large and small accelerators. 



8. Comments on the results 

Others have worked on this subject and there is an overlap between the contents of this 
paper and their work. These previous papers (2 to 6 ) give results for the stopband width 
and for the growth rate inside the stopband. 

The results in this paper include the following. The result given for the tune v near 
the edge of the stopband, z/g, \u — q/2\ = [2\Au\\ue — z/qI]^^^- The result given for the 
beta function near the edge of the stopband, [{j3 — /3o)//3o]max = [2|Az^|/|z^e ~ '^ol]^^^- The 
proof given showing that the real part of v is constant over the stopband at g/2 does not 
depend on perturbation theory, and the result follows from the symplectic properties. The 
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result that all the results found in this paper will also hold for a small accelerator where 
the large accelerator approximation is not used. The result given for the solutions of the 
equations of motion when i^q is inside the stopband, and the proof that the eigenfunctions 
and eigenvalues of the transfer matrix are real inside the stopband. 
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